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Abstract 

Starting from the quantized version of Maxwell's equations for the electro- 
magnetic field in an arbitrary linear Kramers-Kronig dielectric, spontaneous 
decay of the excited state of a two-level atom embedded in a dispersive and 
absorbing medium is studied and the decay rate is calculated. The calcu- 
lations are performed for both the (Clausius-Mosotti) virtual cavity model 
and the (Glauber-Lewenstein) real cavity model. It is shown that owing to 
nonradiative decay associated with absorption the rate of spontaneous decay 
sensitively depends on the cavity radius when the atomic transition frequency 
approaches an absorption band of the medium. Only when the effect of ab- 
sorption is fully disregarded, then the familiar local-field correction factors 

are recovered. 
42.50.-p,42.50.Ct,42.50.Lc 
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I. INTRODUCTION 



Spontaneous emission is a prime example of the action of vacuum fluctuations on phys- 
ically measurable processes. Since the early work of Einstein jTJ spontaneous emission has 
been a major ingredient in the understanding of the effects of what one calls the vacuum in 
quantum field theory. The radiation properties of an excited atom located in free space have 
been a subject of many studies (for a comprehensive list of original articles, see, e.g., ||). 
In particular, the rate of spontaneous emission in free space (half the Einstein coefficient) is 
given by 

_ ^j/i 2 m 

where uj a is the transition frequency of the atom and /x is the dipole matrix element of the 
transition. The question has been arisen of how a surrounding medium modifies that decay. 
Simple arguments based on the change of the mode density suggest that the spontaneous 
emission rate inside a non-absorbing dielectric should be modified according to || 

Tse = riT , (2) 

where n is the real refractive index of the medium. In Eq. @) it is assumed that the local 
field the atom interacts with is identical with the electromagnetic field in the continuous 
medium. Since in reality the atom is in a small region of free space, the local field felt by 
the atom is different from the field in the continuous medium [|J, and the decay rate may 
be expected to be modified to 



SE 



<r , (3) 



where £ is the local-field correction factor. Different models have been used to calculate it. 
In the (Clausius-Mosotti) virtual cavity model it is given by || 
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whereas the (Glauber-Lewenstein) real cavity model leads to 



£cm= , (4) 
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Recently, experiments have been reported |7]|| from which the real-cavity model may be 
favored. 

As already mentioned, in Eqs. (Q) - (|5|) it is assumed that the refractive index of the 
medium, which may vary with frequency [i.e., n — > n(u A ) in Eqs. (1.2) - (1.5)], is real. 
However, in reality the refractive index must be a complex function of frequency, 



n{uj) = T}{u) + i k(uj). 



(6) 



It is well known that causality requires the permittivity of the medium, e(uS) = n 2 (uj), to be 
a complex function of frequency whose real part (responsible for dispersion) and imaginary 
part (responsible for absorption) are related to each other by the Kramers-Kronig relation. 
Only when the atomic transition frequency uja is sufficiently far from a medium resonance, 
so that absorption may be disregarded, the imaginary part of the refractive index (at the 
atomic transition frequency) may be neglected: n(u A ) 

Describing the (undisturbed, continuous) medium in terms of a complex permittivity, in 
PJHI it is argued that Eqs. (|]) - can be extended to the spontaneous emission of an 



atom embedded in a lossy dielectric as 



(7) 



where the local-field correction factors (^) and (^) are now regarded as being squares of 
absolute values, 

n 2 {uo A ) +2 2 
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Further, in 0] the total decay rate is decomposed as 



(9) 



r = r x + 



(10) 



where the rates T 1 - and F% respectively, are related to the transverse and longitudinal 
electromagnetic fields in the medium. The rate T 1 - is identified with the cavity-radius- 
independent rate Tse given by Eq. (0), and it is argued that the rate r", which depends on 
the cavity radius R as r" ~ -R -3 , is responsible for nonradiative decay via energy transfer 
between the atom and the surrounding (absorbing) dielectric. 

^From the study of resonant energy transfer between two guest molecules in a perfect 



lattice of absorbing molecules [II], in [I2| it is argued that (within the approximations 
made) the rate of spontaneous emission is given by Eq. ([?[) together with Eq. @, i.e., with 
the local-field correction factor that corresponds to the virtual-cavity model. However, the 
total decay rate is purely transverse; i.e., it results only from the transverse part of the 
electromagnetic field in the medium, 

r = r- L = rW + r^. (11) 

It consists of an /^-independent far-field term which has the form of Eq. (|7|) [together 
with Eq. (§)] and is interpreted as the spontaneous emission rate, and a /^-dependent term 
T^ 2 \ which in the near-field zone is proportional to R~ 3 and describes nonradiative energy 
transfer. 

Recently it has been shown that the decay rates suggested in P,[lO| for the virtual- 



cavity model are wrong in general, because the quantum vacuum in the presence of a dis- 
persive and absorbing dielectric is not introduced correctly. The fluctuating part of the 
polarization field is not fully included in the local field coupled to the atom and therefore 
effects such as nonradiative energy transfer from the guest atom to the medium via virtual 
photon exchange (i.e., transverse-field-assisted energy transfer) are omitted. It is just the 
contribution to the local field of the fluctuating part of the polarization which gives rise 
to the relevant terms ~ R~ 3 and ~ R' 1 in the transverse decay rate of an excited atom 
surrounded by an absorbing medium IKI. It is worth noting that the results have been 



confirmed within a microscopic approach to the problem more recently 1TJ 



In the virtual-cavity model, the electromagnetic field inside the cavity, i.e., the local 



field, is modified by the presence of the cavity, but the modification of the field outside the 
cavity is disregarded. Hence the local field introduced in this way is not exactly the field 
that couples to the atom in reality. On the contrary, in the real-cavity model the mutual 
modification of the fields outside and inside the cavity are taken into account in a consistent 
way; i.e., the atom interacts with a field that exactly satisfies both Maxwell's equations 
and the fundamental commutation rules of quantum electrodynamics. It may be therefore 
expected that the real cavity model is more suited for describing the spontaneous decay 
than the virtual cavity model. In particular, the Power-Zienau-Woolley transformation (see, 
e.g., ||15|| ) suggests that (in dipole approximation) only the transverse electromagnetic field 
contributes to the decay rate via radiative decay and nonradiative decay associated with 
virtual photon exchange, the latter being typical for an absorbing medium. 

In this article we consider, within the frame of rigorous quantization of the electromag- 
netic field in an arbitrary linear Kramers-Kronig consistent dielectric []i6|-|18|j , the sponta- 



neous decay of an excited atom embedded in an absorbing dielectric, applying the real-cavity 



concept. We find that the rate formulas suggested in JTU] for the real-cavity model are es- 
sentially wrong. At first, only the transverse electromagnetic field contributes to the decay 
rate, i.e., r" =0, which contradicts [ID|. At second, the (purely transverse) rate not only 



contains an /^-independent term but also terms proportional to and -R -3 which are 
closely related to nonradiative decay - a result which also contradicts flnj. As expected, 
nonradiative decay is only observed for an absorbing medium. It is worth noting that when 
the atomic transition frequency is sufficiently far from an absorption band of the medium, 
so that absorption may be neglected, our result exactly agrees with that derived in || for a 
non-absorbing medium. 

The paper is organized as follows. After introducing the quantization scheme, in Sect. |I| 
the problem of spontaneous decay of an exited atom in an absorbing medium is considered. 
In Sect. [TTI| the results for decay rate with the virtual cavity model are outlined, and Sect. [TV| 
presents a detailed analysis of the decay rate with the real cavity model. The results are 
discussed in Sect. [V]. Lengthy calculations are given in the Appendix. 



II. BASIC EQUATIONS 



Our analysis of the spontaneous decay of an excited atom embedded in an absorb- 
ing medium is based on the scheme for quantization of the electromagnetic field in lin- 
ear Kramers- Kronig dielectrics developed in [ — I~B | . We start with the phenomenological 



Maxwell's equations in the (temporal) Fourier space, without external sources, 

V • B(r, uj) = 0, (12) 

V- eo e(r,cu)E(r,cu)] =p(r,w), (13) 

V x E(r, uj) = iuB(r, uj), (14) 

V x B(r, uj) = -i-ze(r, w)E(r, uj) + /i j(r, uj), (15) 

where e(r, uj) = e/j(r, u;) + iei(r,uj) is the (spatially varying) permittivity satisfying the 
Kramers-Kronig relations. When there are no external charges and currents, then p{r,uj) 
and j(r, uj), respectively, are the operator noise charge and current densities that are asso- 
ciated with absorption according to the dissipation-fluctuation theorem. They satisfy the 
equation of continuity, 

V • j{r,uj) = iujp(r,uj), (16) 

- N 

and they are related to the noise polarization P (r, uj) as 

" N 

j{r,uj) = -iujP {r,uj), (17) 
p(r,w) = -V.P (r,uj). (18) 

Let f (r, uj) be an infinite set of bosonic field operators which may be viewed as being collective 
excitations of the electromagnetic field, the medium polarization, and the reservoir. All 
operators in the theory can then be expressed in terms of these basic field operators using 
the relation 
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Ihe 



j(r, u) = u]l — ej(r, u) f (r, w). 



(19) 



In particular, from Maxwell's equations the electric field (in Fourier space) is given by a 
convolution with the classical dyadic Green function, 



E k (r,uj) =i// / d r' ujG kk i(r,r\uj)i k ,(r\u) 



where G kk i(r, r ',u) satisfies the partial differential equation 



dld r k -5 ik [ A r + ^e(r,Lo) 



G kk <(r,r',u) = 5 ik i5{r - r'). 



Integration with respect to to then yields the operator of the electric field as 



E(r) = E (+) (r) + E ( - ) (r), E^( 



E(+) ( 



(20) 



(21) 



(22) 



du E k (r,u) = ifjLo / du d 3 r' uG kk/ (r, r',u)j ,,(r',uj). (23) 
o Jo J - K 



Substituting in Eq. (E3|) for the current density the expression given in Eq. (|T9|) yields the 



electric field in terms of the bosonic basic fields. It can be proven ]I8| that the quantization 
scheme is fully consistent with QED for arbitrary linear dielectrics, i.e., 



4(r),^(r')l = -i%e klm d r J(r-r'), 



(24) 



BJr),B,(r') 



0. 



(25) 



The electric and magnetic fields can be of course expressed in terms of vector (A) and 
scalar (<p) potentials. In what follows we will set the scalar potential equal to zero. This 
gauge condition implies that both the transverse and the longitudinal electric fields are 
obtained from the vector potential 



A(r) = A (+) (r) + A ( -)(r), 



(26) 



4 +) (r) = J^dw J d 3 r'G kk ,(ry,co)i k ,(r',cu). 



(27) 
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Let us now consider the case when an external (two-level) atomic system at position r A 
is present. Treating the interaction of such a guest atom with the electromagnetic field in 
dipole and rotating wave approximations, the Hamiltonian of the total system can be given 
by 

r roc a a 2 ,s r „ ~ -i 

H= d 3 r dwfkufUr,u;)f(r,u) + V fku a A aa - iu 21 A 21 A {+ \r A )-d 21 + H.c. . (28) 
J J o t^x L J 

Here the atomic operators A aa i = \a) (a'\ are introduced, with |a) being the energy eigenstates 

of the guest atom (a = 1, 2). The energies of the two states are ftuj\ and %uj 2 {hu 2 > hu>i), and 

u 2 i = lo 2 — uj\ and d2i, respectively, are the atomic transition frequency and dipole moment. 

Note that in the interaction term in Eq. ( p8|) the A 2 term and the counter-rotating terms 

have been dropped. 

In the Heisenberg picture the equations of motion then read as, on recalling Eqs. ([19]) 
and d?D, 

A 22 = -^ii 21 A(+)(r A )-d 21 +H.c, (29) 
A 11 = -A 22 , (30) 



A 21 = iu 2l A 2l + (r A ) ■ d 21 (A 22 - A u ) , (31) 



/ 4 (r, W ) = -iw/ i (r, W )+ ^^-^^{d 21 ) k GUv A ,v,u)A l2 . (32) 

Substituting in the vector potential in Eqs. (p9|) - (|3l|) for fi(r,u,t) the formal solution of 
Eq. (|32D, i.e., 

(33) 



iJ +) (r,t) = 4 + e L(r,t)+ ^(d 21 ) k [°°du \lmG lk (r,r A ,uj) f dr A 12 (r) 

7rer>cr Jo L Jf 



ireoc* Jo 

a system of integro-differential equations for the atomic quantities is obtained. [Note that 
Eq. ( |A3| ) has been used for deriving Eq. (^3).] At this stage a Markov approximation can 
be introduced, and the integro-differential equations reduce to Langevin-type differential 
equations (Appendix |A|) 

8 



A22 — — r^22 — 

A n = -A 22 , 



A21 ^Ai r + 2(r A ,t)-d 2 i + H.c. 



A 



21 



z(w 21 - Su) - |r] A21 + ^A^(rA,t)-d 2 i (i 22 - i n 



where T is the rate of spontaneous decay of the excited state of the guest atom, 



heoc 2 



Im G kk ,(r A ,r A ,uj A ) 



(34) 
(35) 

(36) 
(37) 



[fi k = (d2i) k , = ^21], an d 5u is the (contribution of the dielectric to the) Lamb shift [see 
Eq. flSfSD I. Note that A£^(r,t) evolves freely. From Eq. (|(| together with Eqs. © and 
it can be proved that the quantization scheme exactly yields, in agreement with the 



dissipation- fluctuation theorem, the relation [19 



Im Gkk> (r, r', a/) 5(u - uj') 



hu 2 



(0| 



|0>. 



(38) 



As long as the Markov approximation applies, the spontaneous decay can be described in 
terms of the rate fl3"7p, the rate formula being valid for arbitrary dielectrics and geometries. 
Especially, for an atom in vacuum we have 



lmG kk/ (r A ,r A ,uj A ) 



6irc 



[see Eqs. (|BT 



5|) for e = 1], which leads to the well-known result ([!]), 

r = r 



(39) 



(40) 



3Ttfieoc 3 

A guest atom in a dielectric is situated in a small free-space region and is surrounded by 
medium atoms. Frequently a cavity model is used for describing the situation. An atom in 
an empty cavity in an otherwise continuous medium is considered and it is assumed that 
the linear dimensions of the cavity are much less then the atomic transition wavelength. In 
particular, for isotropic systems a spherical cavity of radius R may be considered. With 
regard to Eq. (0), the "only" problem that remains is the calculation of (the imaginary 
part of) the classical Green tensor for a dielectric medium of given permittivity which is 
disturbed by a small free-space inhomogeneity. 



III. VIRTUAL CAVITY MODEL 



In the virtual cavity model it is assumed that the field outside the sphere is not modified 
by the small region of free space inside the sphere, and the (local) electric field E'(r, ui) 
inside the sphere is given by [EC 



E(r,cu) = E(r,o;) + — P(r,w) 
3e 



(41) 



where E(r, u) and P(r, ui), respectively, are the electric and polarization fields in the unper- 
turbed continuous medium. From Maxwell's equations (O) and ([15]) together with Eqs. (|i~6|) 



- (19) it is seen that 



P(r, u) = e [e(r, tu) - 1] E(r, u) + P N (r, uo 



(42) 



where 



lhe 



P>,w)=zJ-^ej(r,w)f(r,<<;) 



(43) 



is the noise polarization associated with absorption. For classical optical fields at room 
temperatures the noise polarization weakly contributes to the polarization and the local 
field, and therefore it may be neglected. Obviously, for quantum fields and especially for 
the quantum vacuum, whose coupling to the guest atom gives rise to the spontaneous decay, 
the noise polarization must not be omitted, because it is nothing other but a part of the 
quantum vacuum. Combining Eqs. (|41"|) and ([4^) yields the local-field operator 



M(r, u) = he(r,u) + 2]E(r,u;) + -Lp"(r,w). 
6 6e 



(44) 



It can be shown that the local electromagnetic field satisfies the equal-time commutation 
relations TT3I 



^(r), BI(t')] = -ihe klm d r m 5(r - r') {l + |[e(r, 0) - 1 



(45) 



E' k {vlE[{v') = B> k (r),Bl(r> 



(46) 
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Comparing with the correct commutation relations, we see that the virtual cavity model may 
be regarded as being consistent with QED (over the whole frequency domain), provided that 



e(r,0) < 10; 



(47) 



i.e., the value of the static permittivity must not be too large. 

Now we can turn to the calculation of the spontaneous decay rate, Eq. (|37|) . Recalling 
Eq. (|38|), we may write 

„2 



(48) 



lmG k k'(r,r',uj A ) S(u - u A ) = ( Q \ E k (r,u), E_l(r' ,u A ) |0) 



with \r-r A \,\r' -r A \<R and E' from Eq. (p[). Since E in Eq. (|44]) is determined by Eq. (p0|) 
with the Green tensor for the field in the undisturbed continuous medium, knowledge of the 
imaginary part of that Green tensor is sufficient to calculate the decay rate. However, for 
r, r' — > i\4 a singular contribution to the rate is observed, which reflects the fact that the 
description of the dielectric as a continuous medium contradicts a precise determination of 
the position of the guest atom. The problem might be overcome by regularization, e.g., 
by averaging Eq. ( |Hj| ) over the sphere. Combining Eqs. (|37]) and (|4"8"D and using Eq. (| 
[together with Eqs. (0), (^), and (@)] yields g3| 



CM 



heod 2 



e{uJ A ) + 2 



4 ^A / \ 

+ TjT — ^i\u A )n k n k > Re 

OiI€qC 



e(u A ) + 2 



Gf k ,(r,r',u; A ) 



+ 



9he l 



(the bar introduces averaging over the sphere), where G^(r, r', u;) is the Green tensor of the 
mean field in the undisturbed medium, and e^yi) = e(r A ,u A ). Note that the permittivity 
can be assumed to be constant over the small sphere. The first term in Eq. (|49| ) corresponds 
to the result obtained in RfTCH, without taking account of the contribution of the noise 
polarization to the quantum vacuum. The noise polarization gives rise to the second term 
and the third term in Eq. ( |4TS| ) - terms that are proportional to ei(u A ) and typically observed 
for absorbing media. 
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When the position of the guest atom in the medium is sufficiently far from inhomo- 
geneities (such as the surface of the dielectric body) the Green tensor Gj^.,(r, r', u) in 
Eq. fl49|) may be identified with that for bulk material as given in Appendix [FJ. Insert- 
ing for Gj^,/(r, r', u A ) in Eq. (f49|) the result of Eqs. ([Bl|), (|E2|), and flB5|) and averaging with 
respect to r and r' separately over a sphere, on assuming equidistribution, we derive |2TJ 



CM 



CM 



+ r 



CM) 



(50) 



where Tq M and Tq M , respectively, are related to the longitudinal and transverse parts of the 
Green tensor, 



CM 



r 



4e/(o;^) / c 



27\e(u A )\ 2 \uj a R 



(51) 



e{uo A ) + 2 



f cm = ^o\vM 



+ e I (u A ) [e R (uj A ) + 2] 



15 \uo A RJ 9 



'-k(u a ) 



+ 



2be I (u) A ) 
54 



^) 3 } + °< fl ) < 52 > 



(\Rje(uj A )uj A /c\ <C 1), with T being the free-space spontaneous emission rate defined in 
Eq. ([[]). From inspection of Eqs. (|50| ) - (|52| ) it is seen that, when absorption can be dis- 
regarded, i.e., e T (uj A ) « and hence e(au) ~ e R (uj A ), n{u A ) « ^e R (uj A ), then r CM « 
reduces to Tse given in Eq. (§) with the local-field correction factor (f|). It is further seen 
that for absorbing media the rate Tq M becomes quite different from that given in Eq. (|7|) 
with the local- field correction factor (§), because of the effect of the noise polarization. For 
more details, the reader is referred to ||13|| . Most recently, a more microscopic derivation of 



the decay rate has yielded, apart from regularization factors, the same results ]14 



It is worth noting that the independent terms in Eq. ( |52"D solely result from the noise 
polarization. In particular, the term ~ R~ 3 may be regarded as describing nonradiative 
decay via dipole-dipole energy transfer from the guest atom to the surrounding medium. 
From Eqs. ( p0|) - (^) it is seen that the terms ~ R~ 3 can be combined to obtain an overall 
rate for the nonradiative dipole-dipole energy transfer. Obviously, the decomposition of Tom 
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in Tq M and Tq M has nothing to do with a decomposition in radiative and nonradiative decay 
channels in general. 

It should be pointed out that the averages in Eq. (^9]), which correspond to regularization 
at r— >r', can be taken in different ways. In other words, the independent terms in Eqs. fl5T|) 
and fl52|) are determined only up to some regularization factors. Hence, not only the the 
cavity radius R but also the scaling factors of the absorption-assisted ~ R" 1 and ~ R~ 3 
terms are undetermined in the model. 



IV. REAL CAVITY MODEL 

In the real cavity model the exact Green tensor for the system disturbed by a small 
free-space inhomogeneity is inserted in the rate formula in other words, the electro- 



magnetic field inside and outside the cavity exactly solves Maxwell's equations ( |T2l) - ([15]) 
together with the standard boundary conditions at the surface of the cavity. In contrast to 
the virtual cavity approach, in the real cavity approach the field inside the cavity exactly 
satisfies the fundamental QED equal-time commutation relations ( f24"| ) and (p5|), and the 
Green tensor does not lead to a singular contribution to the decay rate. The Green tensor 
for an inhomogeneous problem of that type can always be written as a sum of the Green 
tensor for a homogeneous problem and some tensor that obeys a source-free wave equation 
and ensures the boundary conditions to be satisfied p2| . Since the guest atom is situated 



in an empty cavity, the relevant Green tensor reads as 

G kk > (r, r a, uj a ) = G^ k , (r, r A , u A ) + G kk > (r, r A , u A ) (r -> r A ) (53) 

where Gj k ,(r, r A , 0Ja) is simply the vacuum Green tensor, which is given by Eqs. ( |B1| ) - 
( p3|) with e{uj) = 1, and G kk /(r, r A , lo a ) describes the effect of reflection at the cavity surface. 
Obviously, G^ k ,(r, r A , u A ) has no longitudinal imaginary part, 

Im G w ( r > r A, u A ) =0 (r -»• r A ) (54) 
Since the tensor G kk >(r, r A , u A ) is related to a source-free problem, it is transverse, and hence 
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Gl kl (r,r A ,u A ) = 



(55) 



The imaginary part of Gkk'{ r A, r A? ma) is therefore equal to the imaginary part of the trans- 
verse part of the Green tensor, so that the rate formula ( |3"?D in the real cavity model reads 



GL 



heoc 2 



Im G kk ,(r A ,r A ,u A ) 



(56) 



In other words, in the real cavity model the longitudinal field does not contribute to the 
decay rate. Thus, the longitudinal decay rate Tq L given in |L0| is an artifact. 

In order to calculate Tgl further, let us again consider a spherical cavity of radius R in 
bulk material, with the guest atom being situated at the center of the sphere. The Green 
tensor for a spherical two-layer system is given in Appendix |C[ From Eqs. (|53| ) - ( |55|) 
together with Eq. (|B5|) [for e(u) = 1] and Eq. ( P22|) it follows that 

u A 



l™-G kk ,{r A ,r Al u A ) 



Qttc 



1 + ReCf (u A ) 



5kk' 



(57) 



with the reflection coefficient C^(uo A ) being given by Eq. ( |C23|) . Hence, for a spherical 
cavity the spontaneous decay rate (|56|) takes the form of 



r GL = r 1 + ReCf^) 



(58) 



where T Q is the free-space spontaneous emission rate ([!]). The reflection coefficient Ci{u A ) 
in Eq. (58) is a function of R and given in Eq. ( p23 ) explicitly. For u A R/c = 2nR/X A <C 1 
we expand it in powers of R to obtain 



3i[e(u A )-l] ( c 



2e(u A ) + l \u A R 
-1 + 0(R), 



9i[4e 2 (cj A )-3e(cu A )-l] / c 



5[2e(u A ) + If 



9e 5 /\u A ) 
uo A Rj [2e(u A ) + If 



(59) 



from which it follows that 



GL 



[ 9e x (u; A ) 




|j2e(a; A ) + l 


2 





9r)(u> A ) 



+ 



9ej(u; A )[28 


e{u A )\ 2 + 12e R {u A ) + l] 


5\2e(uj A ) + 1 


4 



u A R 



|2e(a; A ) + l| 4 
9ft(a; A )ej(a; A ) 
|2e(^ A ) + l| 4 



4|e(c^ A )| 4 + Ae R (cu A )\e(uJ A )\ 2 + e 2 R (cu A ) - ej(uj A ) 
4\e(LU A )\ 2 + 2e R (uj A )} ) + 0(R). 



(60) 
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Needless to say that when setting e(u) = 1, then the free-space spontaneous emission 
rate is recovered. When the atomic transition frequency is far from an absorption band 
of the medium, then absorption may be disregarded, i.e., €i(ua) ~ [and hence €(oja) ~ 



^r^a), ~ \]£r{wa)\- From inspection of Eq. (|60| ) we see that for ei^uja) the 

term proportional to R° is the leading term, which exactly gives rise to the rate formula (|3]) 
together with the correction factor (||), i.e., we recover the familiar result derived in || for 
real refractive index. We further see that for an absorbing medium the rate formula cannot 
be given in the form of Eq. (0) together with Eq. (f|), as is suggested in [pi]]. Equation 
(|60|) reveals that for an absorbing medium terms proportional to R~ 3 and are observed, 
so that the decay rate sensitively depends on the radius of the sphere. In particular, the 
near-field term proportional to R~ 3 can again be regarded as corresponding to nonradiative 
decay via dipole-dipole energy transfer from the guest atom to the medium. 

It should be pointed out that the condition that ujaR/c<^ 1; i.e., the (optical) wavelength 
Xa of the atomic transition must be large compared with the radius R of the cavity, is in 
full agreement with the Markov approximation used in order to introduce a decay rate. 
From inspection of Eq. (|C23| ) it is seen that the (real part of the) reflection coefficient 
(u) becomes a rapidly varying function of frequency for uR/c>l, and hence the Markov 
approximation fails. In that case the sphere acts like a micro-cavity resonator and memory 
effects must be included in the temporal evolution of the atom, which prevents the excited 
state from decaying exponentially. 



To illustrate the results, we have computed the (virtual cavity model) decay rate Tcm, 
Eq. (|(]) - (p2|), and the (real cavity model) decay rate Tq L , Eq. (j60|) , of an atom in a spherical 
cavity of radius R in a surrounding medium with the single-resonance model permittivity 




V. DISCUSSION 



e(u) = 1 + 




(61) 
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Plots of the rates as functions of the atomic transition frequency are given in Figs. [I] - |6|. 
The figures reveal that the two models can yield decay rates that are quite different from 
each other. Far from the absorption band of the medium the difference is rather quantitative 
than qualitative [Figs. and In the absorption band and in the vicinity of the absorption 
band, i.e., in the region between the medium resonance ujt and the longitudinal frequency 
ujl = \JuJt + Wp (in the figures, ujl = 1.1 u>t), a quantitatively and qualitatively different 
behavior of the two rates can be observed [Figs. [[], ||, [5], and ||. In particular, the rate 
obtained with the real cavity model can substantially exceed the rate obtained with the 
virtual cavity model. The differences between the two rates are less pronounced for strong 
absorption; i.e., when the value of the bandwidth parameter 7 in Eq. (|61|) is sufficiently large 
(compare Fig. [I] with Fig. ^, and Fig. [| with Fig. ^). In that region the rates sensitively 
respond to a change of the radius of the cavity (compare Fig. [I] with Fig. ||, and Fig. |3] with 
Fig. d. 

Obviously, an excited atom in an absorbing medium undergoes both radiative and non- 
radiative damping, and in dense media nonradiative decay can be much faster than radiative 
one. In particular, for small cavity radius the ~ R~ 3 dipole-dipole energy transfer terms in 
the two rates can strongly enhance them. Since the radiationless decay typically happens at 
the longitudinal frequency u>l, one observes, for sufficiently small values of 7, a shift of the 
maximum of the decay rate from u>t to lvl with decreasing value of R (compare Fig. |5] with 
Fig. |l[). Even when the atomic transition frequency is relatively far from the medium reso- 
nance, so that the imaginary part of the permittivity becomes relatively small, the values of 
the two rates can notably differ from those obtained from Eq. (0) together with either Eq. (f|) 
or (|^), because of the ~ R~ 3 near field contributions to the rates. It should be stressed that 
Eqs. (^|) - (||) apply only when nonradiative decay can be fully excluded from consideration. 
Otherwise the near-field terms can give rise to observable effects, as is illustrated in Figs. |2| 
and f§ 

The rates r CM and r GL differ essentially in the way the cavity radius is introduced. 
As already mentioned, in the virtual cavity model the needed coincidence limit of the two 
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spatial arguments of the imaginary part of the Green tensor cannot be performed, because 
of the singularity of the Green tensor of the (undisturbed) medium, and regularization is 
required. In the paper, a small fictitious distance |r — r'| 7^ between two neighboring atomic 
positions inside a sphere of radius R is kept in order to get a finite value, and the result 
is then averaged with regard to r and r' separately over the sphere. In contrast, in the 
real cavity model the limit r, r' — > ya can be performed exactly and a proper rate can be 
obtained, R being the radius of the real cavity. ^From the above it is suggested that the 
value of the parameter R may be different in the two models in order to fit each other (note 
that in Figs. H [| the two rates are compared for equal values of R). 

Another consequence of the via smoothing introduced radius of the sphere in the virtual 
cavity model is that there is a non-vanishing ~ R~ 3 longitudinal-field contribution to the 
decay rate. Hence, the nonradiative dipole-dipole energy transfer from the atom to the 
surrounding medium is obtained from the interaction of both transverse and longitudinal 
electromagnetic field components with the atom, T C m = Tcm + ^cm- O n the contrary, 
the real cavity model leads to a decay rate that solely results from the interaction of the 
atom with the transverse field, T GL = Tq L . Here, the dipole-dipole energy transfer fully 
corresponds to a second-order process via virtual photons. It is worth noting that for not too 
small values of the radius of the virtual cavity (in our example, R^O.l Xa) the contribution 
of rjL M to r CM is small, so that it may be disregarded and hence r CM « (see Figs. || 
and H). 

Equation defines the total energy relaxation rate of the (two-level) atom, which re- 
sults from both radiative and nonradiative decay, and the question arises of what is the spon- 



taneous emission rate. In |T(J the transverse contribution to the decay rate is associated with 
spontaneous emission, whereas the longitudinal contribution is associated with nonradiative 
decay. However, the exact result obtained with the real cavity model reveals that there is 
no longitudinal contribution to the decay rate, and hence the transverse contribution must 
be associated with both spontaneous emission and nonradiative decay. Similarly, the decay 
rate obtained from the study of the resonant energy transfer between two guest molecules 
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surrounded by a perfect lattice of absorbing molecules contains only transverse-field con- 
tributions and describes both radiative and nonradiative relaxation processes fliyi^|| . In 
fT2|| it is suggested that the spontaneous emission rate be identified with the i?-independent 
(far-field) contribution to the decay rate. Since the ~ R~ 3 near field contribution may be 
regarded as resulting from nonradiative decay via dipole-dipole energy transfer, the question 
remains of what is the meaning of the remaining terms. Moreover, from our analysis of, e.g., 
the real cavity model it is seen that R must not substantially exceed the atomic transition 
wavelength Xa- Otherwise, the Markov approximation does not apply and the calculated 
decay rate becomes unphysical. In order to answer the question of what is really sponta- 
neous emission, the model should be extended such that light detection at certain distances 
from the guest atom is included. 

Both in the virtual cavity model and the real cavity model the dielectric is described in 
terms of a continuous polarization field that does not resolve the positions of the microscopic 
constituents of the medium. In reality an excited guest atom does of course not interact with 
a continuous medium, but it "sees" the discrete distribution of the microscopic constituents 
of the medium, at least the nearest-neighbor grouping. Hence a refined treatment of the 
medium should also allow for the presence in the cavity of nearest-neighboring medium 
species whose interaction with the guest atom is considered separately. The enlarged cavity 
can then be chosen such that the guest atom cannot "resolve" the discrete structure of the 



medium outside the cavity and the continuous description applies [23,24 
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APPENDIX A: MARKOV APPROXIMATION 



Equation (|3^) can be formally integrated to obtain 



f i (T,u,t) = f beei (T,u,t)+ ^J e -^^(d 21 ) k GUr A ,r,uj) f dr i 12 (r) , (Al) 



where ff r ee(r, u,t) evolves freely. Substituting in the vector potential in Eqs. - fl3T|) for 
fi(r, u, t) the expression given in Eq. (^T) yields a system of integro-differential equations for 



the atomic quantities, which cannot be solved analytically in general. Usually the Markov 
approximation is introduced. It is assumed that (after performing the u integration) the 
time integral effectively runs over a small correlation time interval r c . As long as we require 
that t — t'^>r c , we may extend the lower limit of the r integral in Eq. ( |Al| ) to minus infinity 
with little error. Further we require that r c be small on a time scale on which the atomic 
system is changed owing to the coupling to the electromagnetic field. In this case in the r 
integral in Eq. (|A1|) the slowly varying atomic quantity Ai 2 {r)e %u)2lT can be taken at time t 
and put in front of the integral, 



dr e' iul{t - T) A 12 {t) « / dr e~ iu{t - T) A l2 {r) 



t r t 

a i 12 (t) f dr e -(— )(*-) = A 12 (t) ((u 21 - u) 



(A2) 



[C(x) =ir5(x) +iVx~ 1 ; V denotes the principal value]. Thus, the future of the system is now 
determined by the present time only. We substitute in Eq. ( |A1| ) for the time integral the 
expression given in Eq. (|A2|), calculate the vector potential, Eqs. ( ffijp and (P7|). With the 
help of the relation (see, e.g., jlSlO 

d 3 s — e/(s, u)G km (r, s, u)G* m (r', s, u) = lmG M (r, r', u) (A3) 

we find after some calculation 

A^\r A ,t) = A^ ei (r A ,t) + -^(tfei)* /°°^ ((uo 21 -uj)ImG lk (r A , r A , u) A 12 (t). (A4) 

Txe^& Jo 
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In order to obtain Eqs. (|34] ) - (|36|), we eventually substitute in Eqs. (^) - ([H]) for the 
positive and negative frequency parts of the vector potential the expressions according to 
Eq. (|A4|) . It can be easily seen that the real part of the ( function (i.e., the 5 function) in 
Eq. ([A3]) leads to T given in Eq. (0). The principal- value integral in Eq. ( |X4| ) which arises 
from the imaginary part of the ( function contributes to the Lamb shift and reads 



8ui 



2uj2i(d 2 i)k(d 2 i)k' f°°, lm.G kk/ (r A ,r A ,uj) 
duo 



he c 2 7i 



which can be rewritten as 

2^1i(^2l)fc(^2l)fc' 



8u 



he Q c 2 



ReG k k'(rA,r A ,uj2i) 



UJ - UJ 2 \ 



1 kn G ik {T A , T A , u) 

IT Jo UJ + 0J 2 l 



(A5) 



(A6) 



Equation flA"€| ) holds because of the Kramers- Kronig relation (or Titchmarsh's theorem) for 
the Green function. Note that the real part of the vacuum Green function is infinite for r = 
r' = and regularization is required. The resulting vacuum Lamb shift may be thought of 
as being included in the atomic transition frequency, so that 5uj in Eq. ( |36| ) may be regarded 
as being solely due to the surrounding dielectric. 



APPENDIX B: GREEN TENSOR FOR A HOMOGENEOUS DIELECTRIC 



Following [10,19], the Green tensor for bulk material can be given by 



G kk > (r, r', u) = Gl k , (r, r', u) + G kk , (r, r', u) , 



where (p = r — r') 

G L'( r > r V) 

and 



— S(p) 5 kk > + \ 8 kk i — EhEL j 
3 V P J P 



p 2 J p 3 



' l_ 



:b2) 



kp (kp) 2 (kp) 3 



Ofefc' 



1 3i 
+ 



PkPk' 



kp (kp) 2 (kp) 3 J p 2 



(B3) 
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are related to the longitudinal and transverse electric fields. In Eq. (p3|), the complex wave 
number 



k = Jd^j) - = [ V (u) + m{u)\ - (B4) 

has been introduced. In particular for small values of \kp\, \kp\ <C 1, the exponential e lkp in 
Eq. ( p3|) can be expanded to obtain 

G&(r, v'M = i- 1^ + ^ + ^ fo(u,) + «( W )] <wj + (B5) 



APPENDIX C: GREEN TENSOR FOR AN EMPTY SPHERE SURROUNDED 

BY A HOMOGENEOUS DIELECTRIC 



Following [p5], the Green tensor of a system that consists of an empty sphere surrounded 
by a homogeneous dielectric can be given in terms of spherical Bessel functions and spherical 
harmonics. When r and r' lie in the sphere (with the center of the sphere being the origin 
of the coordinate system), then the associated Green tensor G(r, r', uj) is given by 

G(r, r', cu) = G v (r, r', u) + G(r, r', u), (CI) 

where G v (r, r', c<j) is the vacuum Green tensor, and 

G(r,r>) = ^EEE( 2 '' + \M; 

4ttc £iD n=lm=Q U(n+1) (n+m)\ 



x (2-5 0m ) 



Cf HM :nm ( r, M >m (V, + (u,)N >m (r, N >m (V, ^ 



(C2) 



Here Me nm (r, fc) and N= nm (r, fc) are the (even and odd) vector Debye potentials, and the 
quantities Cff' N (ou) are the generalized reflection coefficients. Introducing the abbreviating 
notations 

J ni = MhR), (C3) 

H ni = hV>{k i R), (C4) 
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1 d[pj n (p)} 



p dp 
I d[ P hW(p)] 
P 



p=k t R 



dp 



p=kiR 



{k\ = \Jt{u) u/c, k 2 — uj/c), the reflection coefficients can be given by 

rpH,V, \r>H,V, \ 
T F,n {u)Rp] n M 



a 



M,N , 



UJ) 



iST(«) 



where 



R v p» 



k 2 H' n2 H n i — k\H' nl H n2 



k2J n iH n2 
k 2 H n2 H' nl 



k\J' nl H n2 
kiH nl H' n2 



k 2 J' nl H n2 - 
k 2 [J n2 H' n2 



k\J n \H' n2 

- J'n2Hn2\ 



F,n(u>) 
V 



k 2 J n iH' n2 - 

k2[J' n2 H n 2 



k\J' nX H n2 

- Jn2H' n2 ) 



^lJ'n2^-n\ ~ 

k 2 [J' n2 H n2 



k\J n2 H' nl 



k 2 J n \H n2 ~ 

k 2 [J n 2H' n2 



k\J n \H n2 

~ J'n2^n2\ 



k2Jn2H' n i — kiJ' n2 H n i 



The vector Debye potentials are defined by 



M. 



(r,k) = V x ipe nm (r,k)r 



Ne nm ( r ,k) = -V x V x [V>m(r,&)r 



with 



/ cos \ 

^>(r, fc) = j n (kr)P™(cos 0) [ ^ j m0, 

and can be given by 

i i \ . /cos\ . . ,, . dP™(cos( 

Me nm (r,k) = —- 3n (kr)P™(cos9)[ . m0e e - j ra (A;r) 

sin 6* \ sin / dU 



+ 



m (r,fc) = ^^ Jn (kr)P;r(cos9)^2)rn<j ) e r 
1 d[r7' n (Ax)] 



At dr 



dP™( cos 9) /cos\ , 

— . m0e eT 

at' \r~ 



sin , 



im 



sin 6 



cos ( 



sin 

.cos. 



22 



j n {kr) is the spherical Bessel function of the first kind and P™(cos#) is the associated 
Legendre polynomial. Note that from Eqs. ( |C14| ) and ( pl5| ) it follows that G(r,r',u) is a 
(two-sided) transverse tensor function. 
Since for fcr->0we have 



3n{kr\ 



kr^O 



(kr) 



1 



(2n + l)!! V 2(2n + 3) 
from inspection of Eqs. (|C171) and (|C18|) we find that 



+ . . . 



(C19) 



M e „ m (r,fc) fc -^°(A;rr, 



(C20) 
(C21) 



Hence, at the center of the sphere only the TM-wave vector Debye potentials Ne 10 (r, k) and 
Ne n (r, k) contribute to G(r,r',cu) in Eq. (|C2|) , 



GkkAr,r',uj] 



r=r'=0 Q n c 



C l ((jj)dkk 1 , 



(C22) 



where [n = 
with 



[i + z(n + 1) — iz 2 n — z 3 n 2 / (n + 1)] e 



sin z — z(cos z + in sin z) + iz 2 n cos z — z 3 (cos z — in sin z)n 2 / (n 2 — V 



(C23) 



Ruj 



(C24) 
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FIGURES 
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FIG. 1. The spontaneous decay rate T, Eq. (|37|), is shown as a function of the atomic transition 
frequency uja near a medium resonance for the model permittivity (|6ll ) {ujp = 0.46 wt, 7 = 0.05u;;r) 
and i? = 0.02Ayi. The solid line corresponds to the real-cavity model, Tgl from Eq. (|60|), and the 
dotted line corresponds to the virtual-cavity model, Tcm from Eq. (|50|), the broken line indicating 



the transverse-field assisted rate Tq M from Eq. (52). 
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FIG. 2. The spontaneous decay rate T, Eq. (|37|), is shown as a function of the atomic transition 
frequency uja far from a medium resonance for the model permittivity (ojp = 0A6ujt, 7 = 

0.05 i^t) and -R = 0.02Aa- The solid line corresponds to the real-cavity model, Tgl from Eq. (|60|), 
and the dotted line corresponds to the virtual-cavity model, Tcm from Eq. (|50|), the broken line 
indicating the transverse-field assisted rate r^ M from Eq. (^). For comparison, the rate Tgl as 
obtained from Eq. (g) together with (|5|) is shown (dashed line). 
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FIG. 3. The spontaneous decay rate T, Eq. (p7\), is shown as a function of the atomic transition 
frequency oja near a medium resonance for the model permittivity ( |6l| ) (top = 0.46cjt, 7 = 0. 2ujt) 
and -R = 0.02A 7 4. The solid line corresponds to the real-cavity model, Tgl from Eq. (|60|), and the 
dotted line corresponds to the virtual-cavity model, Tcm from Eq. (|50|), the broken line indicating 



the transverse-field assisted rate from Eq. (|52j ) 
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FIG. 4. The spontaneous decay rate T, Eq. fl37|), is shown as a function of the atomic transition 
frequency loa far from a medium resonance for the model permittivity (up = 0.46 wt, 7 = 
0.2 ujt) and -R = 0.02Aa- The solid line corresponds to the real-cavity model, Tgl from Eq. (If 



and the dotted line corresponds to the virtual-cavity model, Tcm from Eq. (50), the broken line 
indicating the transverse-field assisted rate Tq M from Eq. (^). For comparison, the rate Tgl as 
obtained from Eq. (||) together with (|5|) is shown (dashed line). 
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ua/d t 

FIG. 5. The spontaneous decay rate T, Eq. (p7|), is shown as a function of the atomic transition 



frequency u>a near a medium resonance for the model permittivity ( plj ) (up = 0.46 wt, 7 = 0.05wt) 
and _R = 0.2Ayi. The solid line corresponds to the real-cavity model, Tql from Eq. (|60|), and the 
dotted line corresponds to the virtual-cavity model, Tcm from Eq. (|50"|), the broken line indicating 



the transverse-field assisted rate Tq M from Eq. (52) 
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FIG. 6. The spontaneous decay rate T, Eq. fl37|), is shown as a function of the atomic transition 
frequency u>a near a medium resonance for the model permittivity (|6l|) (wp = 0.46o;t, 7 = 0. 2wt) 



and i? = 0.2A J 4. The solid line corresponds to the real-cavity model, Tgl from Eq. (60), and the 
dotted line corresponds to the virtual-cavity model, Tcm from Eq. (|50|), the broken line indicating 
the transverse-field assisted rate Tq M from Eq. (|5^). 
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